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ABSTRACT 

Gravitational lensing deflects light. A single lens deflector can only shear images, 
but cannot induce rotations. Multiple lens planes can induce rotations. Such rota- 
tions can be observed in quadruply imaged sources, and can be used to distinguish 
between two proposed solutions of the flux anomaly problem: substructures in lens- 
ing galaxies vs large scale structure. We predict the expected amount of rotation due 
to large scale structure in strong lensing systems, and show how this effect can be 
measured using -mas VLB! astrometry of quadruple lenses with extended source 
structures. The magnitude of rotation is around one degree. The biggest theoretical 
uncertainty is the power spectrum of dark matter on very small scales. This proce- 
dure can potentially be turned around to measure the dark matter power spectrum on 
very small scales. We list the predicted RMS rotation angles for several quadruple 
lenses with known lens and source redshifts. 

Key words: gravitational lensing - cosmology: theory - dark matter - galaxies: struc- 
ture, evolution 



1 INTRODUCTION 

Gravitational lensing results from the deflection of light under the gravitational influence of all 
matter, luminous or otherwise. Its physics is clean, and this eff"ect has allowed the measurement of 
the distribution of dark matter from galaxy-scales (using strong lensing) to the large-scale structure 
of the universe (using weak lensing). 

Usually, many approximations are made to simplify the calculations. One of these is the Bom 
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approximation, where one calculates a small deflection along the unperturbed light path. Some 
effects, such as image rotation due to multi-plane lensing, are not accessible in this approximation. 
Authors have obtained different results for the magnitude of multi-plane weak lensing rotation 
(Jain et al 2000, Cooray & Hu 2002, hereafter CH; Hirata & Seljak 2003, hereafter HS). Schnei- 
der (1997) showed that a strong lens plus spatially constant weak lens system is mathematically 
identical to some other single lens plane system, so only differential rotation at the image positions 
is observable. 

In this paper, we apply the multiple plane lensing calculation to real physical systems: quadru- 
ply imaged quasars. We show in this physical example how the rotation effect can be measured, 
why it is physical and real, estimate its magnitude, and show how it can be used to resolve the 
substructure controversy. 

The rotation of images is of current interests in the context of using gravitational lensing to 
detect substructures predicted by the Cold Dark Matter (CDM) structure formation model. From 
both semi- analytical studies and numerical simulations, it became clear that hundreds of subhaloes 
(substructures) are predicted to exist in a Milky-Way type halo (e.g., Kauffmann et al. 1993; Klypin 
et al. 1999; Moore et al. 1999; Ghigna et al. 2000). In general, about 5-10% of the mass is predicted 
to be in substructures, with a typical mass spectrum of n{M)dM ~ M~^-^dM. If all the substruc- 
tures form stars, then the predicted number of satellite galaxies exceeds the observed number in a 
Milky-Way type halo by a large factor. It is now, however, clear that the correspondence between 
substructures and visible satellite galaxies is not simple (Gao et al. 2004a,b; see also Springel et 
al. 2001; Diemand, Moore & Stadel 2004; Nagi & Kravtsov 2005). In particular, if only some 
substructures house satellite galaxies, then the discrepancy can be alleviated (e.g., Kravtsov et al. 
2004). At present, it is not entirely clear whether the internal kinematics of satellite galaxies are 
consistent with observations (Stoehr et al. 2002; Kazantzidis et al. 2004). Furthermore, the spatial 
distribution of satellite galaxies in the Milky Way is also somewhat puzzling (Kroupa et al. 2004, 
but see Kang et al. 2005; Liebeskind et al. 2005; Zentner et al. 2005). 

One possible way to detect the (dark) substructures is through the gravitational lensing effect. 
Simple analytical models in gravitational lenses often fail to reproduce the observed flux ratios 
(e.g., Kochanek 1991). This discrepancy is commonly referred to as the "anomalous flux ratio 
problem." This has been proposed as evidence for substructures in the primary lensing galaxies 
(e.g., Mao & Schneider 1998; Metcalf & Zhao 2002; Dalai & Kochanek 2002). However, as most 
of the predicted substructures are in the outer part of the lensing galaxies while the lensed images 
are typically at a projected distance of only a few kpc from the center, it is unclear whether the 
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predicted amount of substructures in lensing galaxies by CDM is sufficient (Kochanek & Dalai 
2004; Mao et al. 2004), so it is important to consider other sources of "substructures." Recently, 
Metcalf (2005) proposed that substructures along the line of sight can equally explain the discrep- 
ancy. A key question naturally arises: how do we know that substructures are from the primary 
lens or from elsewhere along the line of sight? 

In this paper, we examine the rotation of images induced by structures along the line of sight 
in gravitational lenses (see also Chen et al. 2003). We use a novel power- spectrum approach and 
consider the fluctuations of surface densities among difl^erent images (separated by few tenths to 
few arc seconds) and their eff'ects on the magnifications. Throughout this paper, we adopt the "con- 
cordance" ACDM cosmology (e.g., Ostriker & Steinhardt 1995; Spergel et al. 2003 and references 
therein), with a density parameter = 0.3, a cosmological constant Q.f^ = 0.7, a baryon density 
parameter Qb = 0.024/z"^, and the power-spectrum normalization erg = 0.9. We adopt a Hubble 
constant of Hq = 70 km s~^ Mpc"^ 

2 LENSING BY LARGE SCALE STRUCTURE IN THE PRESENCE OF A STRONG 
LENS 

Many multiply-imaged gravitational lenses on galaxy-scales have been observed. At the time of 
writing, roughly 100 such systems are known^ The largest systematic survey, the Cosmic Lens 
All Sky Survey (CLASS) has found 22 new galaxy-scale lenses, approximately one half of which 
are quadruple lenses (Browne et al. 2003; Myers et al. 2003). Most of these have high-resolution 
imaging from 0.1" to mas from MERLIN, HST to VLBI. Several of these lenses are resolved into 
multiple components. The system 0128-1-437 provides a good example (Biggs et al. 2004). Each of 
the image has been resolved into three sub-components with VLBI. Such high-resolution images 
provide an excellent test bed for lensing models. 

In addition to the strong lens system which causes the multiple image splitting, all the matter 
along the line of sight will further deflect the light and contribute to distortions in the image. One 
such effect is the apparent rotation of images. It can be shown that a single plane lens has a sym- 
metric amplification matrix, which shears but does not rotate images. With multiple planes, image 
rotation is possible. Observing it may appear non-trivial, since it would require prior knowledge 
of the unlensed image alignment. We will show below (see §3) how this can actually be measured 
if one has a quadruply imaged source with extended structures. 

' see the CaSTLES database: [http://cfa-www.harvarcl.edi]/castles/| 
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Several geometric configurations can lead to image rotation. A strong lens has shear and con- 
vergence of order unity. With sufficiently accurate alignment, a second strong lens could occur 
along the line of sight. Indeed such an example has already been seen - the JVAS/CLASS lens 
system, B21 14+022 (Augusto et al. 2001; Chae, Mao & Augusto 2001), has two lensing galaxies 
at redshifts 0.3157 and 0.5883 respectively, within 2 arc seconds of quadruple radio sources^. 

But in general, the expected variation in surface density due to dark matter is small. Integrating 
the Limber equation, this leads to large scale structure density variations of order a few percent of 
the critical surface density (e.g., Jain et al 2000). So a random cluster lens will typically only have 
weak lenses in its foreground and background. 

The typical splitting angle of the strong lens is around one arc second. Large scale structure 
density fluctuations on such scales are significantly correlated. Since Schneider (1997) had shown 
that perfectly correlated weak lensing screens do not cause observable rotation, one must compute 
the difl^erential weak lensing shear. 

The cross correlation between two image positions separated by angle is defined as 



where ^^(A^) is the two dimensional Fourier (Bessel) transform of Equation (fTTb which will be 
discussed below. The results are shown for four quadruple lenses with known lens and source red- 
shifts, B1422+231 (Patnaik et al. 1992), MG0414+0534 (Hewitt et al. 1992), B1608+656 (Myers 
et al. 1995), and B2045+265 (Fassnacht et al. 1999). We find a ratio of variances between the 
difference of two images at and the individual variances to be 



If the two images are uncorrelated, the difference will have twice the variances of each individual 
image. Figure [l] shows the correlation function vs. image separation. As can be seen, typically we 
have r > 0.5, so the difference mode has a slightly lower variance than each individual image. 

We assume that the various images pass through different parts of the strong lens with corre- 
spondingly different values of shear and convergence, and also differing weak lens deflectors. We 
consider a 2-plane lens, Li is the strong lens, L2 is a weak lens due to large scale structure, O is 
the observer position and S is the source plane. The geometry is shown in Figure El We consider 
a quadratic potential on each lens plane, (p = ax] + IbxiXj + cx\, where (xi, ^2) are the (angular) 
coordinates in the lens plane. The units for the potential are chosen such that 2k = v^0, and the 



(1) 



^ = 2(1 - r). 



(2) 



^ It is present unclear whether the radio sources B and C are lensed images. 
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Figure 1. The correlation function for convergence (k) for different image pairs induced by the large-scale structure; the correlation function at 
zero lag is the variance in k. The boxes indicate the largest and smallest splittings for four CLASS lenses with known lens and source redshifts. 
Because splitting angles are small, the k by the large-scale structure at different image positions is strongly correlated. The variance from large 
scale structure shear at each image position is the correlation function at zero lag. The covariance for images separated by a lag can be seen 
in the plots. Since the covariance (at typical image separations) is of comparable size of to the auto-variance, which leads to a supression of the 
differential shear mode. 



deflection angle is or = Vv0 due to lens, where we use the same notations as in Schneider, Ehlers, 
& Falco (1992). To distinguish the two lens planes, we will use a prime to denote variables in the 
L2 plane. 

For a general anisotropic lens, the deflection angle is a vector. For our quadratic potential, k 
and 7i,2 are constant on the plane, 

K=]^{a + c), 71 = ^(a - c), 72 = b. (3) 

The quadratic potential is the most general function which leads to constant values of k and 7. 
Their constancy on small scales is inferred from the correlation function shown in Figure [U In 
such a potential, the full gravitational lensing effect is straightforward to compute. One can solve 
the full photon trajectory, which are deflected on the lens planes by the gradient of the potential. 

We now denote D\ to be the angular diameter distance from the observer to the first lens, D12 
the distance from the first to the second lens, and the distance from the second lens to the 
source, and so on. We see from Figure [T] that all the relative deflection angles are very small in 
units of radians (~ 10^^), so we can expand to first order in the deflection angle (Schneider et al. 
1992, chapter 9) 
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Figure 2. Lensing geometry. The deflection angles themselves are not observable, but rather the changes in deflection angle. Various quantities are 
indicated in the Figure and used in eqs. I4I7I . 

n=^^i- Du&ii^O - D2,&2{^2\ ^2 = - ^i(^i). (4) 

where r], ^j, ^2 denote the position vectors in the source plane and the first and second lens planes 
(see Fig. El), and the deflection angles in the first and second planes are given as and a2- 
Following Schneider et al. (1992), we define two reduced deflection angles 

^Is ^ £>2s ^ 

(Xi = T^o'i. ^2 = -jz-(^2, (5) 
With these, eq. Q can be recast in a very simple form using only angles: 

y = Xi - Qri(xi) - Qr2(x2), X2 = Xi -/3i2Qri(xi), (6) 

where y = rjlD,, Xj = ^JD^ X2 = ^2/^2, and/3i2 = DnDJ{D2Di^). 

Lensing shear and magnification can be obtained by studying the change of source position y 
resulting from a change of apparent angular position Xi . In particular, the magnification is given 
by (Schneider et al. 1992) 

A = |^ = J-Ui-U2 + A2U2U1, (7) 
axi 

where J is a unit matrix, Ui = dai/dxi, and U2 = da2/dx2. 

We are considering the combined effects of a strong and a weak lens. The contribution of the 
weak lens U2 is small and we are only interested in its contribution to rotation, so we neglect its 
linear efi'ect. In the linear (i.e. Born approximation) regime, it is also observationally not possible 
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to distinguish between contributions from the strong and weak lens. In the product term, we can 
absorb /3i2 into the definition of Uj- We define the critical density for the weak lensing large scale 
structure to be 

(8) 



^crit ^ -Pis 



Ui = 



(9) 



which is the lensing strength of the large scale structure as seen by an observer at the strong lens 
position. This lends a simple observational and computational interpretation of lensing rotation: all 
distortions visible to an observer at the strong lens position enter linearly into the coupling product 
of strong and weak lens. 

An analogous results holds when the weak lensing plane is in front of the strong lens. In our 
notation, we have 

K + 7i 72 ] ( k' + y'l 72 

,72 K-yi ) {72 k' -7\ , 

where k' = and so 

7272 + (1 - + ri)(l -k' + r'l) 72(1 -K-7i) + 72(1 -k' + y[) 
7^(1 -K + yi) + 72(1 -k' - 7'i) 7272 + (1 - - 7i)(l - k' - 7[) 
l-K + 71 72 + 727\ 
72 - 727\ l-K-yi ) 

In the last approximate equality we used the limit that the weak lensing plane (the primed vari- 
ables) are much smaller than the strong lens plane, but keeping the anti- symmetric piece which is 
relevant for rotations. We define the off"-diagonal antisymmetric component as co = y2y\- 

For strong lens systems, the shear on Li is of order unity, while the large-structure shear (y[, 
y'2) is a few percent, corresponding to a rotation of order a degree. Rotation only results when the 
principal axes of the two amplification matrices are misaligned. We choose the coordinates such 
that 
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(10) 
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'7; 72 ^ 


(11) 
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7 J does not contribute to rotation, so we set it to zero. Their product is a pure antisymmetric matrix 



U2U1 = 



772 
-772 ^ 



(12) 



When added to a unit matrix, this corresponds to a rotation matrix by an angle 772. 

In equation (fTUb . we factor the amplification matrix A = AsR((/i), as a product of a symmetric 
matrix and a pure rotation, where the rotation matrix 
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R 



COS0 sin0 



(13) 



A = 



(14) 



- sin cos (p 
We have 

^ \ - K + yi 72 + CO 

72 -co l-K-Ji J 

where a» indicates the importance of rotation. Then we find 

CO 

tan0 = . (15) 

\ - K 

The non-rotating amplification matrix has the same convergence as A up to 0{a?), but has its 
shear components rotated by 0/2. 

The different images pass through different parts of the strong lens, each with its own values of 
shear y. The rotation angle will thus be different for each image. The actual value of the large scale 
structure shear depends on the non-linear power spectrum of dark matter at small physical scales. 
This is both a function of the primordial power spectrum and its slope, and non-linear gravitational 
physics. The most accurate models seem to be a combination of N-body simulations and heuristic 
models based on stable clustering. Calibrations at larger scales indicate consistency at the better 
than 20% level on length scales larger than about 100 kpc. In the standard models, the contribution 
from smaller scales is not dominant, so one would expect forecasts to be good to a factor of two 
(Huffenberger & Seljak 2003). 

To forecast the expected RMS rotation angle, we compute the expected differential variation 
in the weak lens screen shear, and apply to eqs. (I12I15I) . Mathematically, shear is a polarization 
field, which can be described as a trace-free spin-2 tensor field. In two spatial dimensions, any 
trace-free spin-2 tensor field can be decomposed into two kinematic scalars: the "divergence- 
like" component, also known as "E"-mode which is longitudinal to its Fourier decomposed wave 
vector, and a pseudo-scalar "curl-like" or "B"-mode which is transverse to its wave vector. In weak 
gravitational lensing, the "E"-mode is identical to the convergence field /c, and the "B"-mode is 
zero. To compute the statistics of the shear, it thus suffices to calculate the dimensionless variations 
in projected matter surface density. Its variation is given by the Limber equation. It involves the 
projection of a three dimensional non-linear power spectrum 

AHk,z)^^P(k,z) (16) 
to a two dimensional angular power spectrum /(/ -I- l)Ci/2n (Limber 1954; Kaiser 1992, 1998), 
^^^^'^^Q = j £ AHl/Ax(z),z)w(zfx(z)^dz. (17) 



2n I jy dz 
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For the foreground large scale structure, Zi = 0, Zf = Zi, ^x(z) = x(z) where zi is the redshift of 
the strong lens. For the background large scale structure, Zi = Zi, Zf = Zs, ^x(z) = x(zf) - x(z), 
where Zs is the source redshift. We used the Peacock and Dodds (1996) formulation to obtain the 
non-linear power from the linear transfer function given by Bardeen et al (1986). The comoving 
angular diameter distance is 

dz_ 
H{z) 

where H(z) is the Hubble constant at redshift z: 



Da=x{z) = c f-^ (18) 
Jo 



H{z) = //o[(l + zfi^r^z + 1) - ^Kz{z + 2)yi\ (19) 

For the comoving angular diameter distance we used the fitting formula from Pen (1999). In a 
cosmological context, it is convenient to use comoving angular diameter distances and conformal 
time, where light rays propagate as they do in an empty universe (White & Hu 2000). The lensing 
weight is 

w{z) = ^^r.Ho^8{z){l+z) (20) 
where 

g{z) = — — . (21) 

X(Zf)-x(Zi) 

In terms of our previous variables, for the background weak lenses we have = x(z) - xizi), 
Dis = xizs) - xiz), etc. The lensing weighting factor g corresponds to the distance weighted terms 
in Equation dHl). A similar relation holds for the background lenses. 

Tabled] gives the expected rotation angle for a variety of quadruple lenses in the ACDM cos- 
mology. Several simplifying assumptions were made. We attribute all the rotation to the furthest 
image. The change in strong lens y between images is taken to be 0.3, which is multiplied by the 
large scale structure shear in equation (IT^ . What matters is not the change in the absolute value of 
y, but the change in each component. We took the variance of the large scale structure shear to be 
half of the convergence, which corresponds to the variance in 72 in the principal axis frame of the 
strong lens. In practice, only differences in rotation angles are observable, which depends on the 
alignment angles of the shear at diff"erent image positions. There is a contribution to rotation from 
the structure in the foreground as well as the background of the lens. The variances were simply 
added. The ratio of the foreground to background variance is listed in column 5 in Tabled With 
all these caveats, we expect the expected rms image rotation to be good to about a factor of two, 
which is comparable to the expected errors on the theoretical lensing power spectrum. 
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Table 1. Expected image rotation angle for the most magnified image. We assumed a characteristic change in strong lens y of 0.3 between images, 
and K = 0.5 at each image position. Listed are the lens and source redshifts. the expected rotation angle, and the ratio of variances contributed from 
foreground compared to background large scale structure. k+ is the standard deviation of k and y due to foreground large scale structure, while 
is the contribution of the background. 

3 MEASURING ROTATION 

We only consider the simplest case. The source is made up of three components, P^,P^,P^. We 
define P^ to sit at the origin, P^ to sit at {Pi, P^), and P^ at {P\, P^), We can always choose our 
coordinate system this way. 

The lensed image appears at positions A, B, C, D also with three components each. We again 
define the position of the third component to be the origin, and only consider relative distances. 
The apparent position of relative to is P^ = D^A^, and similarly P^ = D^A' where the 
deflection matrix is the inverse of the amplification matrix A defined in equation © for image 
A; each image has its own deflection matrix. It will be convenient to concatenate the two position 
column vectors P^,P^ into a 2x2 matrix P. We can write similar equations for all images, resulting 
in an apparent 16 equations for 16 unknowns: P each has four unknowns, and each symmetric 
amplification matrix has three unknowns. The equations are linear, and homogeneous: 

D^A = P, DbB = P, DcC = P, DoD = P. (22) 

It is clear that one could multiply each solution by a constant and obtain a solution, so one must 
fix one more parameter. Without loss of generality, one could fix P^ = 1, which just fixes a length 
scale. If none of the amplifications are known, the solutions are clearly degenerate between a small 
image that is strongly magnified and a larger image that is less magnified. A further degeneracy 
occurs because we can multiply each equation on the left by an arbitrary shear matrix. Since we 
do not know the intrinsic location of the source substructure positions, this is indistinguishable 
from a constant shear applied to both the lens and the source. This corresponds to a shear plane 
between strong lens and source. For a shear between observer and strong lens, one can simularly 
symmetrize the deflection matrix by multiplying both the lens and the sources by the inverse shear 
matrix. This is in accordance with Schneider (1997). Thus, we only measure the differences in 
shears at the image positions. 

Then we have 16 equations for 15 unknowns, which is over-determined by one. This allows us 
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to solve for one rotation angle. If one assumes the rotation to be dominated by the most magnified 
image, say A, we simply allow to be non-symmetric. This allows one to solve for the rotation 
angle. 

In general, however, the large scale structure shear has two independent components 71,72- 
To make progress, one can assume them to be Gaussian distributed with standard deviation a. 
Their one point function is independent, Pijuji) = N{yi,cr)N(y2,o-). Equation dT^ expresses 
the off-diagonal components of each deflection matrix in terms of the two large-scale structure 
shear components. If we fix a value of 71, we can solve Equation (l22t for 72, giving us an implicit 
definition of 72(71). Then integrating over all possible values of 71 weighted by the probability 
gives the total likelihood for an assumed cr: 



One can then solve for the maximum likelihood value of cr. 
4 DISCUSSION 

If one wishes to observe this effect, several challenges must be overcome. Sources must have at 
least three localizable components, and the position of each must be measured very precisely, to 
better than 1% of the component separation. Positions are also affected by changes in the lens from 
one component position to the next. The latter effect is expected to be significantly smaller than 
rotation, because it depends on the change in shear on small scales, while the rotation depends on 
the large scale structure shear itself. Basically, the rotation is of order of the large scale structure 
shear. The variation of the large scale structure shear on source substructure scales is smaller than 
the shear itself. 

We quantify this as follows: let us assume that the components of the lens have a separation 
of 10 mas. The model assumed that the shear was constant over the apparent size of the source. 
The change in shear across the lens is given by ^^(A^ = lOmas). From Figured] we see that k 
has a differential variance of around 10 at this separation, corresponding to a percent change 
in lensed length scales. This is of comparable magnitude to the rotation effect, which makes it 
desirable to have more than three components to check. In practice, the actual lensing substructure 
is suppressed by several factors. Figure [T] shows the weak lensing correlation in the absence of 
a strong lens. Lines converge behind the lens, which reduces the corresponding physical scales, 
and thus the total variance. And only the variation in shear principle axes across three sub-image 
positions affects the rotation, which is suppressed by another dimensionless factor. 




(23) 
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It is worth noting that HS proved that the rotation oj is identical to the B-mode of the shear, 
up to a constant. In this calculation, we only considered the local value of shears, while the E-B 
decomposition is non-local. We therefore should consider the 4 entries in the magnification matrix 
to be independent. In principle, if one had a very high number density of background sources, for 
example from reionization (Pen 2004), one could also solve for the strong lensing E-B decom- 
posed map (Pen 2000). Such a decomposition also allows one to distinguish between structures 
in the strong lens plane, and contributions from large scale structure along the line of sight. The 
signal to noise required for such an exercise is difficult to achieve with current technology. Using 
quadruply imaged sources allows one to reduce the pairwise noise, and thus solve for the rotational 
component. 

We also note that our calculation did not use the Born approximation. We modeled each lens 
plane by a quadratic potential as inferred from the k correlation function, and computed the full 
deflection trajectories. Recently, calculations of rotations by large scale structure at diff'erent lens 
planes resulted in different answers, depending on whether calculations were done with light rays 
(HS) or with shears on light bundles (CH). In the CH calculation, ray bundles were propagated on 
unperturbed trajectories, but distortions were accumulated. The distortions are relative deflections 
of neighboring rays, which accumulate rotations through the same effect as discussed in this paper. 
This allows CH to see rotation without considering the explicit deflection of the centers of each 
ray. In terms of a full Taylor expansion of light rays along unperturbed paths, this includes some 
of the effects beyond the Born approximation. HS included all leading order corrections to the 
Bom approximation, and thus differed in their answer. In our case, the angular scales involved are 
very small, and we can assume the rotation to be dominated by the common large scale structure 
shear at different image positions. The CH and HS approaches do not differ in our model, since 
the potential is taken to be quadratic, and all second derivatives are constant. 

Our result differ from that of Schneider (1997) in the nature of the effect. Schneider (1997) 
constructs an equivalent single plane strong lens from a strong lens plus constant weak shear. Our 
calculation circumvents the assumption that the weak lens screen is spatially constant. 

Our model assumes a source plane consisting of three compact components with unknown po- 
sitions. If all three components are quadruply lensed, there are a total of 12 angular positions, i.e., 
we have 24 observables. We subtract 8 degrees of freedom for the global lens model, because we 
do not know the macroscopic lens deflection angle. This leaves 16 constraints. Then using Equa- 
tion (0^ we can solve for the relative positions in the source plane modulo a scaling (due to the 
mass sheet degeneracy), which is three numbers, as well as the values of the amplification matrix 
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at each image position, which is four sets of 3 numbers, plus one rotation, for a total of 16 pa- 
rameters. This suggests that small scale structure in a single lens is observationally distinguishable 
from weak shear at a different redshift. 

We have made several approximations, which will affect the results at some level. We used 
the difference in shears at the two furthest image positions. In the analysis, we then assumed that 
three of the images had no rotation, and only the furthest accounted for the multi-plane rotation. 
Reality is more complex, and all images have differential weak lensing shear, and thus some level 
of differential rotation. Simulations are needed to quantify this simplifying assumption. 

In this paper, we have concentrated on the rotation induced by the weak large-scale structure. 
However, for the cases where one has multiple lenses along the line of sight, the rotation can be 
more significant. For the best-fit model in Chae, Mao & Augusto (2001), we find that the two 
Jacobian matrices are given by 



0.176 -0.093 
-0.258 -0.595 



0.888 -0.174 
-0.148 0.481 



(24) 



for images A and D, corresponding to magnifications -7.78 and 2.49 respectively. Clearly for 
the more highly magnified image A, the Jacobian is highly asymmetric and the rotation is quite 
significant. 



5 CONCLUSIONS 

We have computed the expected rotation from uncorrected foreground and background large scale 
structure in strong lensing systems. We have shown that this effect is in principle observable with 
precise VLBI imaging of quadruply imaged lens systems. The rotation is physical and observable 
in the example of a source plane consisting of three point sources. If observed, it can unambigu- 
ously determine if the flux anomaly problem is caused by substructure on the lens plane, or by 
uncorrelated structures along the line of sight. 

In addition, it demonstrates how one can extract information about small scale dark matter 
structure along the line of sight to a lens. A measurement of rotation would measure the variance 
of K from large scale structure at the smallest scales. This affects the scatter in supemovae lensing 
effects, and may have potential to measure the primordial power on small scales, and possibly a 
tilt or running spectral index. 

We thank Peter Schneider for pointing out a flaw in an earlier version of the paper, T. York and 
K.-H. Chae for helpful discussions, and Pengjie Zhang for the lensing Limber code. 
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